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ABSTRACT 


Measurements  of  a  continuous  hydromagnetic  micro¬ 
pulsation  event,  which  occurred  in  the  pre-dawn 
magnetosphere  on  July  29,  1979,  are  reported  as  observed  by 
wave  and  particle  instruments  onboard  the  P78-2  (SCATHA) 
spacecraft.  Calculations  of  the  Poynting  vector  make  it 
clear  that  the  wave  is  a  traveling  Alfven  wave  guided  by  the 
earth's  magnetic  field  line.  Plasma  densities  are 
calculated  at  L  »  7 .  The  phase  relationship  between  the 
plasma  flux  and  the  electromagnetic  fields  is  investigated. 
A  self-consistent  set  of  particle  density  and  particle 
velocity  data  is  determined,  in  order  to  verify  the 
measurement  of  the  electric  field.  Current  magnetospheric 
models  predicting  Alfven  wave  speeds,  Alfven  wave  periods 
and  plasma  densities  are  supported  by  the  study. 
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I.  INTRODUCTION 


A.  GEOMAGNETIC  MICROPULSATIONS 

Geomagnetic  micropulsations  are  low  frequency 
perturbations  of  the  Earth's  magnetic  field.  They  are 
generally  short-lived  variations  whose  amplitude  is  small 
when  compared  to  the  magnitude  of  the  Earth's  main  field. 
In  contrast  to  long  term  variations  of  the  Earth's  magnetic 
field  which  are  of  internal  origin,  geomagnetic 
micropulsations,  as  well  as  magnetic  storms,  are  thought  to 
be  generated  by  solar  forces.  [Ref.  l] 

The  gamma,  defined  as  10~®  Tesla,  is  commonly  used  when 
measuring  these  extremely  small  magnetic  fields. 
Micropulsations  have  amplitudes  ranging  from  fractions  of 
gammas  to  tens  of  gammas  and  periods  ranging  from 
approximately  0.1  second  to  10  minutes.  Micropulsations  may 
be  divided  into  two  main  classes.  Those  with  regular  or 
continuous  patterns  are  called  Pc's  and  those  with  irregular 
patterns  Pi's.  The  system  used  in  this  paper  for 
categorizing  continuous  pulsations  is  described  in  Table  1. 
[Ref.  1] 

B.  CONTINUOUS  PCS  OSCILLATIONS 

Micropulsations  with  the  longest  periods,  PcS's,  may 
have  amplitudes  near  100  gamma,  but  generally  are  seen  to 
have  amplitudes  of  the  order  of  ten  to  twenty  gamma.  A 
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group  of  PcS's  may  last  from  10  minutes  to  several  hours. 
Pc5's  have  teen  noted  to  be  absent  on  days  of  very  low 
magnetic  activity  and  on  very  disturbed  days.  They 

j  generally  occur  during  the  recovery  periods  of  several  days 

1 

which  follow  moderate  magnetic  disturbances.  [Ref.  1] 

TABLE  1 

i  CONTINUOUS  PULSATION  CLASSIFICATION 

K 

[  _ Type  Period  ( sec . ) _ 


Pel 

0.2  - 

5 

Pc2 

5  - 

10 

Pc3 

10  - 

45 

Pc4 

45  - 

150 

PCS 

150  - 

600 

1 .  Previous  Observations 

Magnetic  field  micropulsations  have  been  observed 
from  ground  stations  since  the  19th  Century.  In  1861 
Stewart  observed  what  he  described  as  a  great  magnetic 
disturbance  at  Kew  Observatory  [Ref.  2].  Since  then 
thousands  of  papers  have  been  published  concerning 
micropulsations.  Ground  station  characteristics  of  PcS's 
were  summarized  by  Barfield  and  McPherron  [Ref.  3].  Ground 
stations  report  PcS's  most  frequently  at  dawn  and  dusk. 
These  data  Indicate  that  it  is  common  for  PCS's  to  span  60 
degrees  of  longitude.  Simultaneous  observations  by  magnetic 
conjugate  ground  stations  suggested  a  generation  mechanism 
located  near  the  equator. 

With  the  advent  of  satellite  exploration  in  the  late 
1960 's  the  understanding  of  long  period  micropulsations  was 

8 


'.vVv'v.aV.;.'/.’. 


furthered  with  data  taken  by  satellite  magnetometers. 
Several  studies  were  made  of  data  taken  by  spacecraft  in 
geosynchronous  orbit.  For  example,  correlations  were 
established  between  micropulsation  observations  in  1967  from 
Applied  Technology  Satellite  1  (ATS-1)  and  ground  station 
data  [Refs.  3  and  4].  This  analysis-  showed  that  the  onset 
of  85  percent  of  the  events  occurred  during  the  expans lu.: 
phase  of  magnetospheric  substorms,  suggesting  that  Pc5 ' s 
could  be  generated  by  storm  enhanced  ring  currents . 

In  the  early  1970 's  plasma  scientists  began  to 
correlate  data  made  available  by  multiple  sensors  on 
spacecraft.  Comparison  of  magnetic  field  data  and  high 
energy  plasma  data,  from  Explorer  45,  with  ground  station 
data  showed  that  Pc5  resonance  regions  existed  in  the 
magnetosphere  at  high  altitudes  [Ref.  5].  It  was  suggested 
that  these  resonance  regions  could  be  perturbed  by  a  drift 
mirror  instability.  A  good  review  of  plasma  instabilities 
is  given  by  Chen  [Ref.  6]. 

ATS-6  provided  an  early  example  of  combined  magnetic 
field  and  plasma  flow  measurements  in  1974  [Ref.  7].  Using 
multiple  sensors  it  was  possible  to  calculate  the  Poynting 
vector  using  measured  magnetic  fields  and  inferred  electric 
fields.  This  case  study  showed  the  existence  of  a  standing 
Alfven  wave. 

Multi-sensor  observations  have  become  more  common  in 
the  past  decade.  The  addition  of  mass  spectrometers  brought 


new  dimension  to  magnetospheric  research.  A  recent  example 
is  a  Pc5  observed  by  Dynamics  Explorer  I  in  July  1982  [Ref. 
8].  Using  a  mass  spectrometer,  electric  field  and  magnetic 
field  measurements,  it  was  possible  to  compare  simultaneous 
measurements  of  the  E  X  B  drift  velocities  and  particle 
derived  plasma  drift  velocities  during  a  Pc5  event.  The 
plasma  was  shown  to  be  rotating  perpendicular  to  the 
magnetic  field.  Complete  analysis  of  the  plasma  data  showed 
that  significant  quantities  of  singly  ionized  helium,  oxygen 
and  nitrogen,  as  well  as  hydrogen,  were  present  in  the 
plasma,  which  was  important  in  matching  plasma  theory  and 
observations . 

GEOS-2  data  from  the  double  probe  electric  field 
experiment  have  been  reported  for  PcS's  observed  during  1979 
[Ref.  9].  Observed  correlations  between  wave  amplitude  and 
plasma  density  data  were  used  to  support  current  hydro- 
magnetic  wave  period  models. 

Multi-satellite  observations  can  be  used  to  address 
the  spatial  structure  of  PcS's.  Combined  magnetic  field 
data  from  the  SCATHA,  GOES  2,  GOES  3  and  GEOS  2  satellites 
were  used  to  examine  the  structure  of  a  field  aligned  Pc5 
observed  in  November  1979  [Ref.  10].  The  study  showed  that 
the  wave  had  an  antisymmetric  standing  structure  along  the 
field  line,  and  that  the  wavelength  was  much  shorter  than 
the  field  line.  These  are  characteristics  that  have  not  yet 
been  explained  by  existing  theories. 
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Analysis  of  these  observations,  leading  to  the 
determination  of  wave  modes,  growth  characteristics  and 
propagation  characteristics,  enhance  the  understanding  of 
basic  plasma  physics  and  those  processes  which  determine  the 
dynamics  and  structure  of  the  magnetosphere. 

2 .  Theory 

H.  Alfven  showed,  in  1942,  that  low  frequency 
electromagnetic  waves  can  propagate  in  fluids  of  high 
conductivity  in  the  presence  of  a  constant  magnetic  field 
[Ref.  11]. 

These  waves  are  difficult  to  observe  on  a  laboratory 
scale.  S.  Lunquist  determined  that  for  a  strong  interaction 
to  occur  between  electromagnetic  and  hydrodynamic  systems 
the  criterion  which  must  be  met  is 

BLo(y/p)l/2  >>  jL  (1.1) 

where  B  is  the  magnetic  field,  L  is  the  dimension  of  the 
conducting  fluid,  o  is  the  conductivity  and  p  is  the  density 
of  the  fluid  [Ref.  12].  Since  the  dimensions  of  the 
magnetosphere  are  large,  this  type  of  interaction  is 
possible. 

In  a  Newtonian  fluid,  of  perfect  conductivity  and 
density  p,  small  perturbations  of  the  magnetic  field,  Bq, 
will  propagate  along  the  magnetic  field  line  with  velocity 
given  by 

*  co/k  =  Bq  /  (yoP)l/2  (12) 
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which  is  called  the  Alfven  velocity  [Ref.  6]. 

In  1952,  Alfven  compared  hydromagnetic  waves  to  the 
vibrations  of  an  elastic  string  with  its  ends  fixed  at 
magnetic  field  conjugate  points  [Ref.  13].  A  good  summary 
of  the  analogy  is  given  by  Sugiura  and  Wilson  [Ref.  14]. 
The  derivation  of  this  analogy  is  described  showing  that  the 
plasma  acts  as  the  mass  of  the  string.  The  calculation  of 
perturbation  modes  of  these  string  vibrations,  as  applied  to 
the  earth's  magnetic  field,  is  presented. 

Hasegawa  and  Chen  separate  the  theory  of  long  period 
micropulsations  into  two  categories  and  set  forth  a  theory 
for  each  type  [Ref.  15].  One  type  can  be  explained  as  a 
resonant  Alfven  wave  excited  at  a  local  field  line  by  a 
Kelvin-Helmholtz  instability  at  the  magnetopause.  The 
second  type  is  interpreted  as  an  excitation  of  a  surface 
eigenmode,  a  Kinetic  Alfven  Wave,  excited  by  an  impulse 
energized  ring  current. 

In  this  paper  a  Pc5  event  is  described  which 
exhibits  some  characteristics  of  a  classical  Alfven  wave. 
Its  source  is  believed  to  be  the  ring  current  which 
generates  a  traveling  wave  along  the  magnetic  field  line. 

C.  COORDINATE  SYSTEMS 

The  coordinate  systems  pertinent  to  this  analysis  of  the 
Pc5  data  are  shown  in  Figure  l.  The  magnetic  field  com¬ 
ponents  (Bx/By,B2)  are  described  in  topographic  coordinates. 
X  is  defined  as  North,  Y  points  East,  and  z  points  down. 


V 


T'S’ 


1 

3 


SOLAR  ECLIPTIC 
SYSTEM 


SATELLITE  SYSTEM 


Figure  1. 


Coordinate  Systems 


toward  the  center  of  the  earth,  in  the  radial  direction. 

The  electric  field  data  are  despun  into  two  components 
(Ejj  and  Ej)  of  the  solar  ecliptic  coordinate  system.  The  X 
direction  in  this  system  is  parallel  to  the  earth-sun  line 
and  points  in  the  sunward  direction.  The  Y  axis  is  perpen¬ 
dicular  to  the  earth-sun  line  and  points  toward  dusk.  The  Z 
axis  is  right-hand  perpendicular  to  the  X-Y  plane. 

The  satellite  coordinate  system  is  a  rotating  coordinate 
system  with  X  pointing  in  the  direction  of  the  satellite 
spin  axis,  which  is  essentially  held  constant  throughout  the 
orbit.  The  spin  axis  is  nearly  anti-parallel  to  Y  (SE) .  The 
Y  and  Z  axis  rotate  through  the  plane  perpendicular  to  the 
spin  axis. 

D.  SCATHA  SATELLITE 

The  U.  S.  Air  Force  P78-2  satellite  was  launched  into  a 
near-geosynchronous  orbit  in  January  1979  as  part  of  a  joint 
NASA/Department  of  Defense  program  to  study  Spacecraft 
Charging  At  High  Altitudes  (SCATHA).  A  schematic  of  the 
SCATHA  vehicle  is  shown  in  Figure  2. 

1 .  SCATHA  Orbit  Parameters 

The  SCATHA  orbit  is  elliptical  with  a  perigee  of  5.5 
Re  and  an  apogee  of  7.7  Re.  The  inclination  of  the  orbit  is 
7.5  degrees  and  the  period  is  23.5  hours.  The  satellite 
spin  period  is  59.5  seconds,  with  a  spin  axis  in  the  orbit 
plane  and  perpendicular  to  the  earth-sun  line.  [Refs.  16  and 
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Figure  3  shows  the  satellite's  radial  distance  from 
the  earth's  center  in  earth  radii,  (1  Re  *  6378  km),  for  the 
24  hours  of  Day  210.  The  Pc5  is  observed  from  4  to  7  hours 
UT  during  a  continuous  decline  in  SCATHA's  radial  distance. 
Figures  4  and  5  show  the  geographic  latitude  and  longitude 
during  the  same  period.  During  the  Pc5  event  the  geographic 
latitude  and  longitude  remain  nearly  constant.  The 
variation  in  SCATHA's  latitude  is  a  result  of  the  7.5* 
inclination  in  the  orbit.  The  variation  in  longitude 
reflects  the  fact  that  the  orbit  is  not  geosynchronous,  but 
elliptical.  The  result  of  the  slighly  inclined  elliptical 
orbit  is  a  continuous  change  in  the  satellite's  geographic 
longitude  as  well  as  it's  velocity  and  radial  distance  from 
the  earth.  SCATHA's  velocity  varies  from  2.55  km/s,  at 
apogee,  to  3.74  km/s,  at  perigee.  The  virtue  of  this 
slightly  elliptical  orbit  is  that  the  satellite  can  sweep  a 
wide  sector  of  the  earth's  magnetosphere  while  introducing 
only  small  errors  in  the  raw  magnetospherlc  data. 

In  the  analysis  of  magnetospherlc  data,  the  satel¬ 
lite's  location  relative  to  the  earth's  magnetic  field  is 
important  in  the  classification  of  observed  phenomena.  A 
first  approximation  of  the  shape  of  the  earth's  magnetic 
field  is  that  of  a  dipole  given  by 


R  =  r  sin^e 


(1.3) 


where  r  is  the  distance  from  the  earth's  center  to  the 
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Figure  3.  Radial  Location  of  SCATHA  on  July  29,  1979 
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Figure  4.  Geographic  Latitude  of  SCATHA  on  July  29,  1979 
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intersection  of  the  dipole  field  line  with  the  equatorial 
plane.  R  and  6  are  the  polar  coordinates  of  the  magnetic 
field  line.  The  parameter  used  to  define  the  location  of  a 
magnetic  field  line  is  the  Mcllwain  L  [Ref.  18].  Using  the 
dipole  approximation  for  the  magnetic  field,  the  L  of  the 
field  line  is  equal  to  the  number  of  earth  radii  at  which 
the  field  line  crosses  the  equatorial  plane.  The  value  of  L 
at  other  latitudes  is  given  by 

L  =  Re  /  cos^Oj^  (1.4) 

where  0jjj  is  the  magnetic  latitude. 

A  view  of  the  magnetic  field  dipole  approximation  is 
shown  in  Figure  6  where  field  lines  for  L  =  5,6,  7,  and  8 
are  plotted.  Scatha's  dipole  magnetic  position  between  0400 
UT  and  0700  UT  is  plotted.  At  0400  UT  SCATHA  was  at  L  = 
8.03  and  descended  to  L  »  6.59  at  0700  UT.  It  is  during 
this  period  that  the  Pc5  was  observed. 

The  dipole  approximation  for  the  earth's  magnetic 
field  is  most  exact  at  dawn  and  dusk.  The  approximation  is 
fairly  inexact  at  noon  (where  the  field  is  compressed)  and 
poor  at  midnight  (where  the  field  is  extremely  elongated). 
The  distortion  of  the  earth's  magnetic  field  is  largely  due 
to  the  solar  wind.  A  more  exact  description  of  the  shape  of 
the  magnetic  field  is  given  by  Hess  [Ref.  19]. 
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SCATHA  POSITION  AT  0700  UT 


SCATHA  POSITION  AT  0400  UT 


Figure  6.  Mcllwain  L  and  Geomagnetic  Latitude  of  the  Pc5 
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2 .  Instruments 

The  spacecraft  is  a  cylinder  approximately  1.75 
meters  in  length  and  diameter.  The  instruments  used  in  this 
study  are  the  NASA  Goddard  Space  Flight  Center  electric 
field  monitor,  the  NASA  Goddard  Space  Flight  Center  mag¬ 
netometer  and  the  University  of  California,  San  Diego  (UCSD) 
charged  particle  experiment.  The  locations  of  these 
instruments  are  shown  in  Figure  2.  [Ref.  17] 

a.  The  NASA  Goddard  electric  field  detector 
consists  of  two  antennas  that  are  50  meters  in  length 
projecting  into  the  spin  plane  which  are  used  as  floating 
probes.  The  antenna  wire  is  composed  of  beryllium  copper. 
Each  antenna  is  insulated  except  for  the  last  20  meters, 
which  is  the  active  part.  Differential  signals  between  the 
antennas  give  the  ambient  electric  fields.  The  minimum 
measurable  field  is  0.05  mV/m.  The  time  resolution  of  the 
electric  field  measurements  is  limited  by  the  satellite  spin 
period  (59.5  seconds). 

SCATHA  data  are  despun  into  the  X  and  z 
components  of  the  solar  ecliptic  (SE)  coordinate  system.  By 
this  method  the  two  components  of  the  measured  electric 
field  (Ex  and  Ex),  which  are  in  the  plane  perpendicular  to 
the  satellite  spin  axis  are  known.  It  will  be  shown  that 
the  Ex  component  of  the  field  data  reflects  the  component  of 
the  Pc5  electric  field  which  is  perpendicular  to  the  earth's 
magnetic  field  when  the  satellite  is  near  local  dawn.  The 
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transformation  from  the  measured  fields  to  the  ambient 
fields  requires  the  addition  of  a  0.3  mV/m  sunward  field  to 
account  for  satellite  motion.  [Refs.  16  and  17] 

b.  The  GSFC  magnetometer  is  a  triaxial  flux  gate 
magnetometer  mounted  on  a  4  meter  boom.  The  resolution  is 
0.3  gamma  with  four  vector  measurements  per  second.  Data 
from  the  GSFC  magnetometer  are  used  in  raw  form  and  despun 
into  the  three  vectors  listed  of  the  topographic  coordinate 
system.  [Refs.  16  and  17] 

c.  The  UCSD  charged  particle  experiment  is  an 
electrostatic  analyzer  that  measures  both  ion  and  electron 
flux.  Particle  data  from  the  electrostatic  analyzers  (ESAs) 
are  used  to  measure  the  ion  and  electron  distribution 
functions  and  deduce  the  particle  density,  temperature  and 
velocity.  The  UCSD  particle  detectors  consist  of  five  ESAs 
arranged  in  two  rotating  pairs  of  electron  and  ion 
detectors,  and  one  fixed  ion  detector.  All  three  of  the 
sensors  of  this  instrument  are  contained  in  a  single  package 
mounted  on  the  forward  end  of  the  satellite.  [Refs.  16  and 
17] 

One  of  the  pairs  of  rotating  detectors  covers 
the  energy  range  from  1  ev  to  81  kev  (the  N/S  detectors). 
The  other  pair  (the  E/W  detectors)  and  the  fixed  detector 
cover  the  1  eV  to  1800  eV  energy  range.  The  analyzers  have 
energy  resolution  of  20%  (6E/E)  and  require  16  seconds  for 
one  64  step  energy  scan.  [Refs.  16  and  17]  A  more  complete 
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II .  OBSERVATIONS 


A.  PARTICLE  AND  FIELD  DATA 

During  the  morning  hours  of  Day  210  SCATHA  recorded  a 
disturbance  in  the  earth's  magnetic  field,  beginning  at  L  = 
7.96,  as  the  spacecraft  approached  the  maximum  inclination 
of  its  orbit,  7.5  degrees.  A  strong  magnetic  field 
micropulsation,  a  Pc5,  was  recorded  on  each  of  SCATHA 's 
sensors . 

Figure  7  is  an  equatorial  plane  plot  of  SCATHA 's  orbit 
on  July  29,  1979  (Day  210).  The  location  of  the  extremes  of 
the  SCATHA  orbit  are  shown  with  respect  to  the  earth  sun 
line.  The  sector  of  the  orbit,  in  which  the  Pc5  event  was 
observed,  is  also  indicated  in  Figure  7. 

1 .  Magnetic  Field  Data 

Figure  8  shows  the  SCATHA  magnetic  field  data 
recorded  from  0300  UT  to  0800  UT.  Included  are  the  total 
magnetic  field,  Bq,  and  the  three  topographic  components  of 
the  field  (Bjj,By,B2).  The  significant  rise  in  the  magnitude 
of  Bq,  the  total  field,  and  Bj^,  the  northerly  component,  is 
due  to  the  continuous  descent  of  the  spacecraft  through  the 
earth's  magnetic  field  for  that  period  of  it's  orbit. 
During  that  time  period,  SCATHA' s  radial  distance  diminishes 
from  7.3  Re  to  5.3  Re,  and  L  decreases  from  8.28  to  6.16. 
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Figure  7.  Equatorial  Plane  Plot  of  SCATHA  Orbit 
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Figure  0.  Magnetic  Field  Data  fror  :300  to  08 


a»4 


EAST-WEST  MAGNETIC  FIELD  (  By  ) 


Figure  9.  The  Pc5  in  the  East-West  Component 

of  the  Magnetic  Field 


The  amplitude  of  the  radial  component  of  the 
magnetic  field,  Bg/  varies  from  l  to  4  gamma  during  this 
period.  The  magnitude  of  this  component  is  largely  due  to 
the  tilt  of  the  earth's  main  magnetic  field  vector  in  the  Z 
direction.  The  perturbations  in  are  a  result  of  a  small 
component  of  the  Pc5  in  that  direction.  The  Pc5  is  observed 
in  the  plot  of  the  east-west  component  of  the  field.  By, 
beginning  at  4:16  UT  and  continuing  until  just  before  dawn, 
06:56  UT.  Since  the  perturbations  in  B^  and  B2  are  small, 
the  Pc5  event  is  largely  contained  in  By. 

The  By  data  are  expanded  in  Figure  9  for  the  time 
period  0400  to  0700  UT  (2:14  to  5:21  local  time).  A  charac¬ 
teristic  of  this  Pc5,  which  makes  it  fairly  unique,  is  the 
packetized  nature  of  the  pulsations.  These  packets  imply 
different  regions  of  hydromagnetic  resonance  in  the  magneto¬ 
sphere.  The  packets  are  labelled  A  through  E  for  con¬ 
venience  in  comparison  to  the  electric  field  data  which 
follows.  The  period  of  the  pulsations  decreases  with  time. 
It  will  be  shown  that  the  period  is  a  function  of  L  and  of 
the  particle  density. 

2 .  Electric  Field  Data 

Figure  10  is  a  plot  of  the  electric  field  data 
measured  by  SCATHA  from  0300  to  0800  UT.  As  mentioned 
earlier,  the  electric  field  data  is  despun  into  two  com¬ 
ponents  of  the  solar  ecliptic  (SE)  coordinate  system.  The 
Ejj  component  points  in  the  sunward  direction  and  remains 


DAY  210  ELECTRIC  FIELDS 


almost  constant  at  1.5  mV/m  throughout  the  event  (2  mv/m  has 
been  added  to  the  value  of  E^,  in  Figure  10,  in  order  to 
separate  from  E2).  This  component  is  a  result  of  the- 
SCATHA's  relative  motion  through  the  plasma,  the  static 
convective  electric  field,  the  sunward  electric  field  caused 
by  photoemission  [Ref.  17]  and  a  partial  projection  of  the 
Pc5  electric  field  perturbation  vector  in  the  X  direction. 

The  E2  component  is  perpendicular  to  the  solar 
ecliptic  plane.  This  component  is  nearly  equal  to  the  total 
electric  field  component  of  the  Pc5.  A  -0.43  mV/m  offset 
seen  in  the  data  is  attributed  to  error  fields  and  a  steady 
state  convective  field. 

Figure  11a  is  an  expanded  plot  of  E2  from  0400  UT  to 
0700  UT.  The  packetized  nature  seen  in  the  magnetic  field 
pulsation  is  reflected  in  the  electric  field.  For  con¬ 
venience,  the  electric  field  data  is  seperated  into  named 
resonance  packets  as  listed  in  Table  2. 

TABLE  2 

HYDROMAGNETIC  RESONANCE  PACKETS 
Packet _ t  begin _ t  end 


(hours 

UT) 

A 

4.41 

4.90 

B 

5.03 

5.46 

C 

5.65 

6.02 

D 

6.10 

6.40 

E 

6.63 

6.95 

The  component  of  the  electric  field,  Ep,  which 
results  from  a  perturbation  in  By  (topographic  coordinates) 
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ELECTRIC  FIELD  (  Ez  ) 


FIG.  11a 
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Flgvire  11.  Measured  Electric  Field  and  Low  Energy  Ion 
Distribution  Function  Showing  Modulation 
of  Ion  Flux  by  the  Field 
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should  be  perpendicular  to  both  By  and  Bq.  (SE 
coordinates)  is  believed  to  be  a  projection  of  Ep  onto  the 
electric  field  antenna.  This  is  based  on  the  assumption 


that  there  can  be  no  component  of  the  electric  field 


sustained  in  the  direction  of  Br,. 


6E  •  Bq  =  0 


(2.1) 


The  electrons  in  a  plasma  are  free  to  move  along  the 


magnetic  field  line  and  quickly  reduce  any  potential 


difference  which  might  be  impressed  along  that  direction, 


Figure  12  shows  the  geometry  for  extracting  Ep  from 
E^.  Here  a  is  the  angle  between  B^  and  Bq  and  also  the 


angle  between  the  radial  direction  and  Eq.  p  is  the  sum  of 


the  satellite's  geographic  latitude  and  the  satellite's 


inclination  from  the  solar  ecliptic  plane.  The  angle  ® 
between  E^  and  Ep  is  given  by 


®  «  90  -  (a  +  p) 


, (2.2) 


Ep  -  E2  /  COS  S 


(2.3) 


Table  3  lists  these  angles  and  the  correction 


factors,  e,  for  calculating  Ep  for  the  hours  of  0400  to  0700 
UT.  The  worst  case  value  is  only  60  percent.  The 


correction  factor  in  Table  3,  e,  is  used  as  follows: 


e  =  1  /  cos  ® 


lEpI  -  lEj 


(2.4) 
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Fig’are  12.  Calculation  of  E. 


TABLE  3 

ELECTRIC  FIELD  DATA  CORRECTION  FACTORS 


t  (UT) 

a 

P 

S 

e 

(degrees ) 

0400 

54.4 

19.2 

16.4 

1.04 

0430 

50.1 

18.0 

21.9 

1.08 

0500 

47.8 

16.7 

25.5 

1.11 

0530 

41.3 

15.2 

33.5 

1.20 

0600 

37.6 

13.4 

39.0 

1.29 

0630 

32.9 

11.5 

45.6 

1.43 

0700 

28.7 

9.2 

52.1 

1.63 

The  factor  which  most  affects  the  magnitude  of  this 
correction  is  the  shape  of  the  earth's  magnetic  field  at 
different  locations  in  the  orbit.  It  can  be  seen  in  Figure 
12  that,  for  small  a,  less  of  the  perpendicular  field  is 
projected  onto  the  antenna.  Near  midnight  the  field  is 
extremely  elongated  giving  a  large  angle  between  Bq  and  Bjj, 


and  thus  a  small  angle  between  E  and  the  Z  antenna.  At  dawn 
the  magnetic  field  most  closely  approaches  a  dipole  and  so 
this  angle  is  small.  During  the  analysis  of  the  electric 
field,  with  this  small  correction  in  mind,  the  actual 
electric  field  data  will  be  assumed  to  represent  the  perpen¬ 
dicular  electric  field  (E2  -  Ep) . 

3 .  Low  Energy  Ions 

Figure  11b  shows  a  plot  of  the  log  (base  10)  of  the 
distribution  function  of  the  39.7  ev  ions  measured  from  4  to 
7  hours  UT.  The  Pc5  can  be  seen  in  the  ion  flux.  A 


comparison  of  this  plot  with  the  electric  field  shows  that 
the  39.7  ev  ion  flux  is  modulated  in  phase  with  the 


perturbation  of  the  electric  field,  6E2.  This  reflects  the 
convective  drift  Induced  in  the  ions  by  the  electric  field. 
Higher  energy  ions  and  electrons  are  also  modulated,  as 
shown  in  the  expanded  presentation  of  the  D  packet  which 
follows . 

B.  A  CLOSER  LOOK,  THE  D  PACKET 

Since  the  D  resonance  packet  exhibits  the  largest 
amplitude  fields,  it  was  chosen  for  closer  analysis.  This 
packet  begins  at  6.1  hours  UT  and  ends  at  6.4  hours  UT.  At 
this  time  SCATHA  is  at  7.3*  North,  having  recently  passed 
the  maximum  inclination  of  the  orbit.  The  D  packet  is 
observed  for  a  period  of  17.23  minutes.  The  center  of  this 
region  is  at  L  =6.97  and  at  4:21  local  time.  During  this 
time  SCATHA 's  average  velocity  was  3.25  km/sec.  After 
taking  into  account  the  decrease  in  SCATHA 's  radial  position 
through  this  period,  this  resonance  packet  appears  to  span 
an  arc  of  3300  kilometers, 

1 .  Field  Relationships 

Figure  13  shows  plots  of  spline  fitted  E2  and  By 
versus  time  for  the  D  packet.  Within  the  resonance  region, 
E  leads  B  by  an  average  of  130*  ±10*,  (-  3tt/4),  and  the 
average  period  of  the  pulsations  is  224  seconds.  Averages 
are  used  since  the  periods  and  phases  are  continuously 
changing,  as  if  the  fields  are  driven  by  a  multi- frequency 
or  non-harmonic  force.  The  electric  field  amplitude  varies 
from  0.62  to  1.8  mV/m.  The  magnetic  field  amplitude  varies 
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from  6  gamma  to  12  gamma.  The  large  amplitudes  at  the 
center  of  the  packet  distinguish  the  center  of  resonance  in 
this  region. 

The  electric  field  vector  and  tht,  magnetic  field 
vector  of  the  Pc5  lie  in  a  plane  perpendicular  to  the  main 
magnetic  field  line.  The  electric  field  vector  pulsates 
toward  and  away  from  the  earth  in  a  direction  between  the 
radial  direction  and  the  (SE)  +Z  direction.  The  magnetic 
field  vector  pulsates  in  a  nearly  east-west  direction.  If 
the  fields  were  in  phase  the  wave  would  be  propagating  in 
the  -X  direction  (south  along  the  magnetic  field  line).  A 
standing  wave  will  exhibit  n/2  phase  difference.  A  phase 
difference  of  -tt  would  indicate  the  wave  is  propagating 
northward . 

The  average  phase  of  fields  in  Figure  13  is  3tt/4. 
This  implies  that  the  wave  is  not  a  pure  propagating  wave. 
The  continuously  changing  phase  of  the  fields  and  the  motion 
of  the  spacecraft  through  different  regions  of  resonance 
make  the  phase  interpretation  somewhat  ambiguous. 

The  complexity  of  the  phase  relationship  is  born  out 
by  a  phase  hodogram  plot  of  versus  6By  in  the  D  packet. 
Figure  14.  In  a  hodogram  a  propagating  wave  would  give  a 
diagonal  line,  and  a  standing  wave  a  circle  (if  E  and  B  are 
plotted  in  Volts/meter  and  Tesla).  This  hodogram  implies 
that  the  wave  is  propagating  at  times,  where  the  plot  is 
nearly  linear.  At  other  times  it  shows  the  phase  of  a 
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standing  wave,  where  the  plot  becomes  more  circular.  The 
average  phase  of  37t/4  is  reflected  in  the  overall  elliptical 
nature  of  the  plot. 

E  and  B  are  related  by  E  *  v^j^B.  Hence,  the  slope  of 
a  line,  or  the  inclination  of  an  ellipse,  in  Figure  14  gives 
The  slope  of  the  major  ellipse  is  0.153  (mV/m) /(gamma) 
which  gives  a  phase  velocity  of  153  km/sec. 

2 .  Particle  and  Field  Relationships 

Particle  count  rates  were  measured  by  the  UCSD 
electrostatic  analyzers  and  the  energy  distribution 
functions  were  calculated.  Ion  and  electron  fluxes  are 
obtained  for  63  energy  levels  every  16  seconds. 

Figure  15  contains  plots  of  the  log^^Q  F  versus  time 
for  low  energy  ions,  high  energy  electrons,  and  two  sets  of 
high  energy  ions  measured  between  6.1  UT  and  6.4  UT.  The 
data  sets  chosen  were  from  energy  levels  which  exhibited  the 
highest  particle  count  rates.  This  was  done  in  order  to 
minimize  the  statistical  error. 

A  plot  of  log  F  versus  time  shows  the  fluctuations 
in  the  energy  of  the  ions  with  time.  The  main  cause  of  the 
change  in  the  kinetic  energy  of  the  low  energy  ions  is  the 
EX  B  force,  which  is  generated  by  the  Pc5  wave.  In  the  D 
packet  region  the  magnetic  field  strength,  Bq,  is  approxi¬ 
mately  140  gamma.  Hence,  a  2.6  mV/m  electric  field  should 
give  the  plasma  an  E  X  B  drift  velocity  equal  to 

Vj5  «  6E2  X  Bq  /  1Bq2  ]  -  18.6  km/sec  (2.5) 
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Here  the  upper  limit  of  the  electric  field  corrected  with 
Table  3  values  is  used  (5E  •  e). 

If  the  plasma  flow  is  perturbed  by  a  change  in  the 
fields,  then  the  kinetic  energy  is  given  by 

U  -  1/2  m  (V  -  vj))^  (2.6) 

-  1/2  m  (V2  -  2  V  Vj)  cos  e  Vj)2) 

where  V  is  the  particle  and  v^  the  E  X  B  flow.  If  v  >>  vp 
(which  is  true  for  40  ev  ions)  and,  if  0  -  0,  then  Equation 
(2.7)  can  be  expanded . 

U  -  1/2  m  (V  -  vj3)2  (2.7) 

-  1/2  mv2  -  (mV  vjj) 

If  the  particle  energy  distributions  were 
Maxwellian,  they  could  be  represented  by 

F  -  n  (m/2TrkT)3/2  exp  -(U/kT)  (2.8) 

and  by  Equation  (2.7) 


log  F  -  constant  +  (log  e)  (mv/kT)  •  vj) 

Hence,  the  plot  of  log  F  for  40  ev  ions  is  propor¬ 
tional  to  vj).  Comparision  of  39.7  ev  ions  with  the  electric 
field  in  Figure  15a  shows  this  to  be  the  case.  The  ion  flux 
is  out  of  phase  with  the  electric  field  by  tt  (-E^  is 
plotted).  The  phase  difference  with  the  electric  field  is 
due  to  the  unidirectional  counting  of  the  ESA's,  since  they 
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are  looking  in  the  direction  of  the  positive  EX  B  force, 
the  count  rate  is  maximum  when  SEg  is  negative.  Thus,  the 
low  energy  ion  flux  is  modulated  by  the  change  in  the  flow 
energy  induced  by  the  E  X  B  force.  The  higher  energy 
particles  respond  primarily  to  B  [Ref.  22]. 

Figure  15b  shows  that  9.3  kev  ions  are  effectively 
modulated  in  phase  with  the  magnetic  field  data.  These 
higher  energy  particles  are,  of  course,  unaffected  by  the 
small  change  in  the  velocity  caused  by  the  E  X  B  force. 
Hence,  the  changes  in  9.3  kev  ion  flux  reflect  the  effective 
motion  of  the  different  magnetic  field  lines  with  respect  to 
the  satellite.  Figure  15c  shows  4.7  kev  electrons  which,  in 
like  manner,  show  flux  modulations  in  phase  with  the 
magnetic  field. 

At  higher  energies  the  modulation  of  the  ion  flux 
begins  to  lag  the  magnetic  field.  This  is  shown  in  Figure 
15d  where  the  flux  of  the  47.4  kev  ions  lags  the  magnetic 
field  by  an  average  of  7r/2.  This  difference  reflects  the 
large  gyroradius  (225  km)  of  these  high  energy  particles. 

These  comparisons  illustrate  the  point  that  the 
plasma  is  moving  with  the  magnetic  field  lines  of  force. 
Hence,  the  plasma  and  the  field  lines  oscillate  together  as 
if  the  particles  are  frozen  to  the  field  lines  [Ref.  6]. 

3 .  Plasma  Parameters 

The  plasma  density,  temperature  and  flow  velocity 
can  be  determined  by  assuming  that  ion  distribution 


functions  are  nearly  Maxwellian  and  fitting  Maxwellian 
distributions  to  the  data. 

Analysis  of  the  ion  data  requires  a  multi-species 
population  due  to  the  various  masses.  The  distribution 
function  becomes 

F  -  S  =  2  nj^  (mj^/2TtkT) exp  -(U^/kTi)  (2.9) 

and  log  F  for  each  population  is 

log  F  =  log  {  n  (m/2TrkT) )  -  (U/kT)  log  e  (2.10) 
.  -  f(p,T)  -  (log  e  /  kT)  U 
-  f(p,T)  +  f(T)  u 

To  begin  the  fitting  process  the  distribution 
function  of  the  ion  data  is  calculated  and  a  plot  of  log  F 
versus  energy  is  made.  As  an  example.  Figure  16  shows  a 
plot  of  the  ion  distribution  function  calculated  from  ion 
flux  recorded  by  the  N/S  detector  at  6:12  UT.  This  time  was 
selected  because  it  is  near  a  large  amplitude  in  the 
electric  field.  It  is  here  that  the  highest  drift 
velocities  can  be  expected  and  the  effect  of  the  positive 
satellite  potential  shadowing  the  lowest  energy  ions  (H"*") 
the  smallest. 

In  this  fitting  process  the  plasma  was  assumed  to 
contain  H'*',  He"^  and  0"^  ions  only.  Once  the  ion  species  are 
assumed,  there  are  three  unknowns  for  each  species  in  the 
distribution  function.  The  temperature,  kT,  the  density,  n. 
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Figure  16.  UCSD  Ion  Distribution  Function 
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and  the  flow  velocity,  v.  A  self  consistent  set  of  values 
is  required.  In  particular,  consistent  flow  velocities  are 
required  since  all  ions -move  together.  One  additional  free 
parameter  is  the  satellite  potential. 

The  temperature,  kT,  is  approximated  to  be  10  ev. 
In  a  more  distinctly  separated  energy  distribution,  each  ion 
species  would  show  a  well  defined  peak  and  the  width  of  the 
peak  would  define  the  temperature.  In  the  present  case,  the 
curvature  of  the  distribution  function  is  a  result  of  both 
temperature  and  velocity. 

Since  the  detector  cannot  be  assumed  to  be  looking 
in  a  direction  exactly  parallel  to  the  flow,  a  look  angle, 
<t),  must  be  included.  The  choice  of  this  angle  affects  the 
calculated  density  and  flow  velocity.  In  this  case  10 
degrees  was  chosen. 

A  peak  in  the  distribution  function  determines  a 
velocity,  Vp,  given  by 

Vp  -  (2U/m)^''2 

-  Vq  cos  <t)  (2.11) 

where  Vp  is  the  plasma  drift  velocity.  Hence,  the  error 
introduced  in  the  calculated  velocity  by  an  erroneous 
selection  of  the  view  angle  is  proportional  to  cos  <J>.  An 
error  of  15*  in  the  approximation  of  (J>,  will  cause  an  error 
in  Vp  of  less  than  10  percent.  Olsen  discusses  this  process 


and  shows  that  the  error,  introduced  in  the  calculation 
of  n  is  more  severe  [Ref.  23]. 

C  «  exp  -(V£)2  -  cos^^j)  (2.12) 

In  Figure  16,  the  peaks  occur  at  approximately  8  ev, 
35  ev  and  100  ev.  The  lower  energy  peak  is  assumed  to  be 
due  to  H'*',  the  middle  energy  peak.  He"'',  and  the  highest 
energy  peak  ,0'*’.  An  approximation,  Vp,  for  the  drift 
velocity  of  each  population  can  be  made  from  the  measured 
energy  peak. 

Since  V  =  (2U/m)l/2  (2.13) 

Vp  39  km/ s 
Vp  jjp  41  km/s 
Vp  0  -35  km/s 

In  the  example.  Figure  16,  the  data  are  fitted  with 
a  Maxwellian  distribution  function  for  H"*",  for  He"''  and  for 
O'*".  Consistency  requires  that  the  Maxwellian  plots  for  all 
three  species  have  the  same  drift  velocity  and  a  reasonable 
density  ratio.  The  values  used  for  this  fit  are  35  km/s 
velocity,  a  10  ev  temperature,  a  look  angle  equal  to  10 
degrees,  and  ion  number  densities  1.0,  0.5,  and  0.3  (cm“3) 
for  H'*’,  He"*"  and  O'*'.  The  results  of  this  fitting  process 
using  UCSD  data  over  the  D  packet  are  listed  in  Table  4. 


TABLE  4 

PLASMA  ION  NUMBER  DENSITY  (1-500  ev) 

Ion _ n _ V _ 

(cm“^)  (km/sec) 


H+ 

1.27 

± 

.67 

35  ± 

5 

He+ 

0.38 

± 

.05 

35  ± 

5 

0+ 

0.30 

± 

.10 

35  ± 

5 

Total 

1.95 

± 

.73 

-  35 

The  flow  velocity  of  35  km/sec  contains  both  the  flow 
induced  by  the  PC5  and  the  steady  background  flow.  It  is 
estimated  that  the  background  convective  flow  is  l  -  10 
km/sec.  The  steady  state  electric  field  responsible  for 
such  a  flow  (0.2  -  2.0  mV/m)  is  poorly  measured  by  the 

electric  field  Instrument.  The  remaining,  pulsation  induced  •. 
flow,  is  approximately  25  km/sec. 

From  Equation  (2.6),  the  SEg  X  Bq  velocity  was 
calculated  to  be  18.6  km/sec.  This  implies  a  steady  state 
flow  velocity  near  16.4  km/sec.  Since  6E  -  Vjj  •  Bq  and  Bo  = 
140  gamma,  a  convective  electric  field  as  large  as  2.3  mV/m 
is  implied  from  the  electric  field  data.  Analysis 

throughout  the  D  packet  is  consistent  with  this  interpre¬ 
tation. 

An  alternative  interpretation  is  that  6E  is  twice 
the  value  measured  by  the  floating  probe.  Experience  with 
the  electric  field  data  by  the  principal  investigator,  T.  L. 
Aggson,  has  suggested  problems  with  the  amplitude  of  the 
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measured  fields.  An  aggreement  within  a  factor  of  two,  at 
this  point,  is  considered  to  be  supportive  of  the  direct 
measurements  of  E. 

The  complete  ion  distribution  is  needed  to  determine 
the  total  ion  density.  The  high  energy  UCSD  data  were 
combined  with  the  high  energy  ion  composition  data  from  the 
SC8  High  Energy  Ion  Composition  Experiment  [Ref.  21].  SC8 
determined  average  ion  densities  in  the  0.1  to  32  kev  range 
for  the  one  hour  period  from  0600  to  0700  UT.  The  average 
H'*’  density  was  (1.974  cm"^},  the  He'*’  density  was 
(0.089  cm“^}  and  the  O'*'  density  was  (0.398  cm”^}.  The  SC8 

data  are  shown  in  Figure  17.  The  line  plot  in  Figure  17  is 
a  plot  of  log  F  for  the  total  population.  The  low  energy 
(100  ev  -  500  eV)  0"^  ions  are  the  same  population  inferred 
by  the  fitting  process  of  the  UCSD  data  in  Figure  16.  The 
UCSD  (0  -  80  kev)  ions  are  shown  in  Figure  18.  The  10  keV 

ions  are  part  of  the  bulk  of  the  plasma  sheet  population, 
while  the  47  keV  ions  are  from  the  tail  of  the  distribution 
function.  The  SC8  data  are  combined  with  the  parameters 
obtained  from  the  UCSD  experiment  and  reported  in  Table  5. 

TABLE  5 

PLASMA  ION  DENSITIES  (.01  -  32  KEV) 


Ion 

n 

( cm" 3 ) 

H+ 

3.24  ±  .68 

He+ 

0.47  ±  .15 

0+ 

0.40  ±  .17 

TOTAL 

4.1  ±  1. 

a  [0+] 


>  CH+] 


ENERGY  (KEV) 
X  [HE  +  ] 


F  (PUVSMA) 
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4 .  Electron  Density 

The  calculation  of  the  ion  density  depends  largely 
on  the  distribution  function  of  the  ion  population  being 
nearly  Maxwellian.  An  error  in  the  look  angle  can  introduce 
a  large  error  in  n.  Since  the  number  density  of  the 
electrons  should  be  equal  to  the  number  density  of  the  ions 
in  a  plasma,  the  value  of  ng  is  a  good  check  of  the  ion 
density  calculations. 

All  the  members  of  the  electron  population  have  the 
same  mass,  and  that  mass  is  known.  Hence,  one  of  the 
unknowns  that  hindered  the  ion  calculation  is  removed.  A 
useful  characteristic  of  magnetospheric  electrons  is  that 
sections  of  the  distribution  function  plot  have  clearly 
defined  slopes  which  identify  different  energy  populations 
of  electrons. 

In  a  plot  of  the  electron  distribution  function 
versus  energy 

log  Fg  =  log  (ng  (mg/2TTkTe)  3/2)  _  log  e  (l/kTg)  Ug 

=  c  +  b  Ug  (2.14) 

where  b  is  the  slope  of  each  section  of  the  distribution 
function  and  c  is  the  y  intercept.  From  Equation  (2.14) 

kTg  =  -log  e  /  h  (2.15) 

ng  -  10^=  /  (mg  /  2tt  kTg)  3/2 
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Figure  19  is  a  plot  of  the  electron  distribution 
calculated  from  UCSD  data  taken  at  6:12  UT.  The  units  of 
Figure  17  are  (sec^/km®)  versus  (keV).  Correcting  for  units 

(2.16) 

n©  =  (10^  X  10"24j  /  (9.05  X  10“^3  /  kT  (ev)  }3/2 

In  Figure  19  two  energy  sets  of  electrons  are  identified. 
By  the  slope  intercept  method,  the  17.5  kev  electrons  are 
shown  to  have  a  number  density  of  0.007  cm“^.  The  number 
density  of  the  1.65  kev  electrons  is  shown  to  be  1.77  cm"^. 
By  expanding  the  lower  energy  section  of  this  data  two  other 
electron  groups  were  identified.  The  number  densities  which 
were  calculated  for  each  set  of  electrons  is  reported  in 
Table  6 . 


TABLE  6 

ELECTRON  NUMBER  DENSITY 


15  - 

80 

17.50 

0.007 

± 

.001 

1.0  - 

15 

1.65 

1.77 

± 

.09 

0.2  - 

1.0 

0.71 

1.46 

± 

.09 

0.0  - 

0.2 

0.12 

0.67 

± 

.04 

TOTAL  3.9  ±  .2 


Hence,  the  calculated  ion  number  density.  (4,1  cm~^) 
is  consistent  with  the  independently  determined  electron 
density. 
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5.  Self  Consistency  of  the  Data 


It  is  now  important  to  determine  whether  these 


numbers  represent  a  self-consistent  relationship  between  n, 
B,  E  and  V.  Since  the  plasma  acts  as  a  parallel  resistive 


circuit  on  the  electric  field  antenna,  a  major  concern  is  to 


verify  the  amplitude  measurements  of  E2. 

Cummings  used  the  following  equation  to  deduce  the 
density  of  the  ambient  plasma  from  field  and  particle  data 
measured  by  the  ATS-6  spacecraft  in  geosynchronous  orbit 


n  -  6b2  /  viq  <ni>  v^^  (2.17) 

where  <m>  is  the  average  particle  mass  [Ref.  7].  This 
relationship  is  used  to  investigate  the  consistency  of  the 
ion  density  calculations. 

Assuming  that  the  plasma  is  made  up  of  He'^  and 
O'*",  and  using  the  values  from  Table  6,  the  ion  mass  density 
is  given  by 

p  ■  S  nj^  mj^  (2.18) 

-  (3.24(1)  +  .47(4)  +  .40(16))  mp 
p  -  (11.52  mp)  kg/cm^ 

and  <m>  -  p  /  n  «  2.81  mp 

mp  ■  1.672  X  10“^^  kg. 


where 


Equation  (2.17)  is  then  rewritten, 

v2  «  68^  /  Mo  ^ 

V  -  6B  /  (Mo 

Since  I V  I  -  1 6E  /  Bq I 

6E  -  6B  Bq  /  (Mo  (2.19) 

-  6.4  X  10“^  (6B  Bq)  mV/m 
where  the  units  of  6B  and  Bq  are  gamma. 

In  the  D  packet  Bq  ■  140  gamma 

6  <  6By  <  12  gamma  (2.20) 

Therefore  the  magnitude  of  the  electic  field  inferred  by 
mass  density  is 

5.4  <  6E  <10.8  mV/m  (2.21) 

Ihe  measured  values  of  6E2  were 

0.62  ^  6E2  i  1.8  mV/m  (2.22) 

Taking  into  account  the  electric  field  correction  factor  of 
1.43,  from  Table  4,  the  magnitude  of  the  perturbed 
electric  field  is  calculated  to  be 

0.9  <  6Ep  <  2.6  mV/m  (2.23) 

These  values  are  a  factor  of  3  lower  than  predicted  by 
Equation  (2.21).  A  factor  of  2  could  be  accounted  for  if 
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the  impedance  of  the  plasma  was  comparable  to  the  internal 
impedance  of  the  electric  field  experiment  (this  is  roughly 
the  case).  A  small  error,  not  larger  than  a  factor  of  2,  in 
the  electric  field  data  is  implied  by  the  comparison. 

6 .  The  Alfven  Velocity 

The  same  comparison  can  be  made  in  terms  of 


velocities.  The  Alfven  velocity  is 


Va  «  w/k  =  Bq  /  (Mo  P) 


(1.2) 


Since 


then 


Curl  E  »  -dB/dt 

ik  Ep  -  -(-iW)  By 
6Ep/By  -  w/k  = 


(2.24) 


Using  upper  and  lower  bounds  of  the  field  amplitudes 
given  by  Equations  (2.20)  and  (2.23),  the  Alfven  velocity 
calculated  from  measured  E  is 


328  km/sec  <  <  378  km/sec 


Using  the  ion  density  relationship 


(2.25) 


'^A  "  ®o  (^^o  “  500  km/sec 


(2.26) 


The  electric  field  data  appears  to  be  low  by  a  factor  of  3 
in  this  comparison. 

7 .  The  Poynting  Vector 

The  nature  of  the  propagation  of  the  PCS  can  be  seen 
by  examining  the  relationship  of  6E2  and  6By.  If  the 


..->  •.'.  .V  L' 


micropulsation  was  a  classical  plane  Alfven  wave,  then 

S  -  (l/y)  E  X  B 

-  (1/yo)  6E  X  6B 

-  1  (l/Mo)  (5Ey  6Bz  -  6E2  6By) 

which  assumes  the  other  components  of  the  Poynting  vector 
are  small.  6Ey  is  not  measured,  but  &B2  is  small,  so 

S  -  1  (l/Po)  (-5Ez  6By)  (2.27) 

In  Figure  20  Px  -  (6E2  6By),  so 

S  -  -i  (I/Mq)  (Px)  *  10~^2  watts/m^  (2.28) 

The  negative  values  of  Px  in  Figure  20  indicate  that 
the  wave  is  propagating  in  the  +X  direction,  northward  along 
the  magnetic  field  line.  The  wave  must  have  been  generated 
south  of  7*  N.  The  fact  that  Px  is,  at  times,  zero  or 
positive  indicates  that  the  micropulsation  is  not  a  simple 
traveling  wave,  but  is  superimposed  on  a  smaller  wave 
traveling  in  the  southward  direction.  This  is  the  mixed 
standing  /  propagating  nature  inferred  previously  in  Figure 
14. 

And  finally,  the  magnitude  of  the  Poynting  vector  is 
calculated 

Px  -  (1.8  mV/m)  (-12  gamma) 

iS|  -  (l/Vo)  |Px|  X  10~^2  (2.29) 

-  17.2  yWatts/m^ 
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The  Poynting  vector 
(Calculated  in  the  D  Packet) 


Figure  20. 


This  presumes  that  the  electric  field  measurement  is 
correct.  Given  the  possible  error  in  6E,  the  magnitude  of 
the  Poynting  vector  could  be  as  high  as  50  pwatts/m^. 

8 .  Dependence  of  the  Period  on  L 

The  period  of  the  Pc5  depends  on  the  length  of  the 
field  line  and  the  Alfven  velocity  of  the  wave  [Ref.  24]. 
The  variation  of  micropulsation  period  with  latitude  can  be 
expressed  as  a  function  of  L  and  n,  the  plasma  number 
density.  [Ref. 25]. 

t  -  (2Tr  Re/®e)  •  (1^0  n/Q)^/2  (2.30) 

Bg  «  3.12  X  10“^  Tesla 

Q  is  proportional  to  an  Alfven  wave  eigenvalue  which  is 
calculated  for  different  values  of  L,  density  gradient,  and 
wave  type.  In  a  simpler  form  this  equation  corresponds  to 

t  •  27T  (m/k)^'^2  , 

The  SCATHA  data  can  be  used  to  determine  the  dependence  of 
the  period  of  this  Pc5  on  L.  If  the  period,  t,  is 
proportional  to 

log  T  -  log  C  +  m  log  L  (2.31) 


where  C  is  a  constant. 

Figure  21  shows  the  average  micropulsation  period  in 
each  resonance  paclcet,  plotted  against  L.  The  slope, 
determined  by  linear  regression,  is  3.96  ±  .07.  Hence,  the 
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PCS  PERIOD  VERSUS  L 


THE  PERIOD  OF  A  PCS  WAVE 
IS  PROPORTIONAL  TO 


LOG  [  L  ] 


Figure  21.  Pc5  Period  Versus  Mcllwain  L 


SCATHA  data  is  in  close  agreement  with  the  theoretical 
relationship 

t  «  (2.32) 

Having  verified  the  L  dependence  of  the  period.  Equation 
(2.30)  is  rewritten  and  utilized  to  predict  a  mass  density 
to  compare  with  the  previous  measurements. 

p  -  n  <m>  =  (X  Be  /  2n  )2  x  (Q  /  Mq)  (2.33) 

For  the  values  x  ■  224  sec  (in  the  D  packet) 

L  =  7.0 

Be  -  3.12  X  10“5  Tesla 

n  <m>  -  2.514  mp  •  Q  (cm"^)  (2.34) 

Orr  and  Matthew  calculate  Q  for  a  toroidal  mode  wave 
in  plasma  density  gradients  which  obey  power  laws  between 
and  r2  [Ref.  25].  At  L  -  7  these  values  are 

3.6  (for  r4)  <  Q  <  4.2  (for  r2)  (2.35) 

The  mass  density  Inferred  by  Equation  (2.33)  is 

9.05  mp  <  n  <m>  <  10.56  mp  (cm"^)  (2.36) 

The  measured  mass  density  was  11.52  mp  (cm'^).  If  the 
theoretical  relationship  is  valid  for  this  Pc5  (it  has  been 
shown  that  it  is  substantively  valid  by  the  agreement  with 
the  power  law  of  the  period),  then  the  difference  between 
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predicted  particle  density  and  measured  density  implies  that 
the  mass  density  has  been  slightly  overestimated.  However, 
a  deviation  of  the  wave  period  of  only  10  seconds  from  the 
average  period,  used  in  Equation  (2.34),  would  bring  the 
predicted  mass  density  into  exact  agreement  with  the 
measured  mass  density. 

Since  the  values  for  D  were  derived  from  a 
theoretical  eigenvalue  equation  for  the  toroidal  wave, 
complete  aggreement  with  the  mass  density  results  of  this 
study  is  not  demanded.  However,  this  comparison  shows 
strong  agreement  with  both  variables,  L  and  p.  The  plasma 
density  gradient  is  the  unknown  variable. 

C.  SOLAR  AND  MAGNETIC  FIELD  CONDITIONS 

Geomagnetic  and  Solar  Data,  published  by  the  Journal  of 
Geophysical  Research,  indicates  that  Day  210  was  a 
moderately  active  day,  Kp  »  3  to  5,  with  no  magnetic 

substorms  recorded.  Forty  eight  hours  prior  to  the  Pc5 
event  there  was  a  period  of  very  high  magnetic  activity  in 
which  a  12  hour  magnetic  storm  occurred  [Ref.  26].  These 
conditions  support  the  statement  by  Jacobs  [Ref.  1]  that 
Pc5 ' s  generally  occur  during  recovery  periods  following 
magnetic  disturbances. 

D .  OTHER  OBSERVERS 

SCATHA  observed  the  Pc5  at  331  degrees  east  longitude  on 
L  -  7.  Magnetograms  from  the  National  Institute  of  Polar 
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Research  from  the  conjugate  ground  stations  at  Leirvogur, 
Iceland  (68.2*  N  338.3*  E)  and  at  Syowa  Base,  Antarctica 
(69.0*  S  39.6*  E)  indicate  that  this  event  was  observed  at 
these  stations  simultaneously  from  0600  to  0700  UT.  The 
time  resolution  of  the  ground  station  data  was  not  clear 
enough  to  say  more.  These  stations  are  located  at  the 
invariant  latitude  for  L  =  7.  The  simultaneous  observation 
of  this  event  would  imply  that  the  generation  was  occuring 
near  the  equator. 

There  is  no ' indication  of  this  observation  at  Narssarss- 
uaq  Station,  Greenland  (61.2*  N  314.2*  E)  or  at  St.  Johns 
Station  (47.6*  N  307.3*  E) .  Magnetograms  from  the  NOAA 
GOES  2  and  GOES  3  satellites  were  also  scanned.  They  showed 
frequent  magnetic  activity,  but  none  which  could  be 
correlated  with  the  Pc5.  These  satellites  were  in 
geosynchronous  orbit  at  252*  East  and  225*  East 
respectively. 

The  lack  of  Pc5  activity  at  locations  west  of  the  SCATHA 
sector  indicates  a  limit  on  the  temporal  and  spatial  extent 
of  the  event.  The  Pc5  appears  to  have  been  confined  to  45 
degrees  of  longitude  in  the  morning  sector  observed  by 


SCATHA. 


III.  DISCUSSION 


A.  THE  CHARACTER  OF  THE  PCS 

In  Chapter  II,  the  observable  quantities  characterizing 
this  PCS  event  were  investigated  by  focusing  on  data  taken 
at  the  temporal  and  spatial  center  of  the  event.  The  wave 
was  shown  to  be  a  traveling  Alfven  wave  guided  by  the 
magnetic  field  line.  The  Alfven  velocity  was  nearly  900 
kilometers  per  second.  A  calculation  of  the  Poynting  vector 
showed  the  wave  was  traveling  northward. 

The  electromagnetic  field  perturbation  was  shown  to  be 
occuring  in  a  plane  perpendicular  the  the  magnetic  field 
line.  The  magnetic  field  lagged  the  electric  field  by  37t/4. 
The  amplitude  of  the  magnetic  field  perturbations  was  near 
12  gamma.  The  measured  electric  field  reached  2.S7 
millivolts  per  meter.  By  requiring  consistency  of  the  data, 
the  measured  electric  field  was  shown  to  be  low  by  a  factor 
of  2  to  3 .  The  period  of  the  PcS,  at  L=7,  was  224  seconds 
and  the  period  of  the  wave  was  shown  to  be  proportional  to 
l4. 

The  plasma  density  was  estimated  to  be  1.9  x  10“^^ 
(gm/cm^).  The  ion  number  density  was  determined  to  be  4.1 
(cm~^).  This  is  consistent  with  previous  magnetospheric 
studies.  The  calculation  of  the  Alfven  velocity  *  900 
km/sec  is  also  consistent  with  previous  studies.  Most 
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recently,  empirical  models  of  average  magnetospheric  plasma 
and  field  data  were  developed  to  study  Alfven  and  Fast  Mode 
wave  speeds  [Ref.  27].  These  models  predict  plasma 
densities  of  3-5  (cm“3),  plasma  temperatures  of  30-100  ev 
and  Alfven  wave  speeds  of  500  -  1000  km/sec  at  L  =  7.0  in 
the  morning  sector  of  the  magnetosphere. 

B.  SPATIAL  EXTENT  OF  THE  PC5 

The  most  interesting  factor  in  this  Pc5  event  is  the 
packetized  nature  of  the  micropulsation  resonances  through 
this  sector  of  the  magnetosphere.  The  breadth  and  location 
in  these  regions  is  presented  in  Table  7 . 


TABLE  7 

RESONANCE  PACKET  LOCATION  AND  SIZE 


t  ID  Width  Latitude  L  Local  T 


4.66 

A 

5300 

7.38 

7.8 

2:45 

0.45 

5.25 

B 

5300 

7.50 

7.5 

3:23 

0.50 

5.84 

C 

4200 

7.45 

7.2 

3:54 

0.32 

6.26 

D 

3300 

7.29 

6.9 

4:21 

1.00 

6.80 

E 

3600 

6.96 

6.7 

4:55 

0.72 

7.46 

F 

5400 

6.11 

6.4 

5:42 

0.49 

The  r  factor,  in  Table  7,  is  a  normalized  resonance 
strength  for  each  region,  which  is  the  ratio  of  the  electric 
field  amplitude  in  a  region  to  the  electric  field  in  the  D 
packet.  These  resonance  regions  are  seperated  by  sectors  of 
low  resonance  which  vary  in  distance  from  1200  kilometers  to 
3600  kilometers. 


C.  UNRESOLVED  ISSUES 


A  complete  understanding  of  this  event  would  require  a 
vector  description  of  the  data.  There  are  significant 
aspects  of  the  Pc5  which  were  neglected  by  this  study.  One 
is  the  effect  that  the  nutation  of  the  spacecraft  spin  axis 
has  on  the  raw  data.  The  nutation  effect  can  be  seen  as  a 
small  pulsation  in  each  component  of  the  magnetic  field  data 
in  Figure  8.  Another  aspect  is  the  significance  of  the 
smaller  component  of  the  Pc5  which  is  in  the  radial 
direction,  B2.  The  component  has  an  average  amplitude  of 
4  gamma.  This  component  has  a  small  effect  on  the  amplitude 
of  the  Pc5  magnetic  field,  but  has  a  significant  effect  on 
the  rotational  motion  of  the  field.  A  study  that  included 
this  component  would  provide  a  more  realistic  description  of 
the  Pc5  event. 

A  most  difficult  aspect,  which  is  left  for  the  student 
who  follows  this  study,  is  the  separation  of  the  spatial 
character  of  the  event  from  the  temporal  description  which 
this  study  provides.  Many  of  the  effects  seen  in  the  data, 
such  as  the  sinusoidal  variation  of  the  resonance  in  each 
packet,  have  been  described  but  not  explained. 

D .  RECOMMENDAT IONS 

1.  The  spatial  and  temporal  extent  of  this  Pc5  could 
not  be  determined  by  SCATHA  alone.  Specifically,  it  is  not 
known  if  local  dawn  was  a  boundary  of  the  Pc5  or  that  the 
Pc5  ended  at  all  locations  at  6:56  UT.  The  radial  boundaries 
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of  the  event  are  also  left  unresolved.  The  observation  of 
the  Pc5  ended  when  SCATHA  passed  L  =  6.6.  This  may  have 
been  a  coincidence,  but  it  could  be  a  major  characteristic 
of  this  phenomenon.  A  correlation  with  a  satellite  that  was 
East  of  SCATHA  could  provide  an  answer  (such  data  may  be 
available  from  the  European  GEOS-2  satellite).  Better  time 
resolution  of  the  conjugate  ground  station  data  would 
determine  the  location  of  the  generating  mechanism. 

2.  Approximate  densities  were  calculated  in  the  region 
of  highest  resonance,  the  D  packet.  The  wealth  of 
information  made  available  by  SCATHA  provides  an  area  for 
further  study  into  the  magnetospheric  plasma  profile  between 
5.5  and  7.5  Re. 

3.  The  theoretical  models  presently  available  describe 
standing  and  compressional  Alfven  waves.  Since  this  event 
provides  extensive  data  from  a  traveling  wave  event,  further 
investigation  into  this  area  of  micropulsation  theory  is 
warranted . 

4.  Theorists  have  recently  shown  great  interest  in 
understanding  the  resonance  conditions  for  Alfven  waves.  As 
emphasized  in  this  paper,  the  packetized  nature  of  the 
resonance  conditions  in  the  morning  sector  of  the  magneto¬ 
sphere  is  displayed  vividly  by  the  SCATHA  data.  These 
observations  should  be  used  to  improve  the  understanding  of 
Alfven  wave  resonance. 
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IV.  CONCLUSIONS 


Observations  of  a  Pc5  event  occuring  between  L  -  6.6  and 
L  -  8.0  in  the  pre-dawn  sector  of  July  29,  1979,  at  7*  N, 
were  reported.  Hydromagnet ic  resonance  regions  spanning 
3000  to  5000  kilometers  with  centers  separated  by  2400 
kilometers  were  identified. 

The  plasma  density  in  the  highest  resonance  region  was 
shown  to  be  1.92  x  10“^^  gm/cm^.  Ion  flux  at  40  ev  was 
shown  to  be  in  phase  with  the  electric  field  .  Ion  fluxes 
at  9  kev,  and  5  kev  electrons,  were  shown  to  be  in  phase 
with  the  magnetic  field.  The  47  kev  ions  lagged  the 
magnetic  field  by  ti/2  radians. 

Particle  data  from  the  UCSD  electrostatic  analyzers  was 
shown  to  be  an  effective  means  of  cross  checking  electric 
field  data.  The  electric  field  data  was  assumed  to  be  a 
projection  of  the  perturbation  electric  field  onto  SCATHA's 
antenna.  By  comparing  this  data  with  ion  density 
predictions,  the  amplitude  of  the  measured  electric  field 
data  was  shown  to  be  a  factor  of  2-3  too  low  to  be 
consistent  with  the  complete  set  of  data  (this  is  considered 
to  be  a  small  error). 

The  periods  of  the  micropulsations,  between  L  =  6.3  and 
7.7,  were  shown  to  be  proportional  to  the  fourth  power  of  L, 
which  supports  present  models  of  resonant  Alfven  wave 


activity  in  the  magnetosphere.  Alfven  velocities  were  shown 
to  be  nearly  900  kilometers  per  second  within  high  resonance 
regions  of  the  event. 

The  electromagnetic  pulsations  were  in  a  plane 
perpendicular  to  the  magnetic  field  line.  In  the  Pc5,  the 
magnetic  field  lagged  the  electric  field  by  3tt/4.  The 
Poyntlng  vector  was  shown  to  be  in  the  direction  of  the  main 
magnetic  field  vector,  confirming  that  the  wave  was  a 
traveling  Alfven  wave.  Simultaneous  observations  by 
conjugate  ground  stations  indicated  that  the  generation 
could  be  occurring  at  the  equator. 

The  unique  packetized  resonance  regions  exhibited  by 
this  PCS  as  well  as  the  non-classical  wave  character  of  this 
event  indicates  that  further  theoretical  study  utilizing 
data  from  SCATHA  is  warranted. 
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